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Abstract 

We perform a systematic study of the discrete time Quantum Walk on one dimension using 
Wigner functions, which are generalized to include the chirality (or coin) degree of freedom. In 
particular, we analyze the evolution of the negative volume in phase space, as a function of time, 
for different initial states. This negativity can be used to quantify the degree of departure of the 
system from a classical state. We also relate this quantity to the entanglement between the coin 
and walker subspaces. 
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I. INTRODUCTION 



Quantum Walks (QW) are considered a as a piece of potential importance in the design of 
quantum algorithms ll- 5 1, as it is the case of classical random walks in traditional computer 

n 

science. As in the case of random walks, QW's can appear both under its discrete-time [6fl or 
continuous-time Q form. Moreover, it has been shown that any quantum algorithm can be 
recast under the form of a QW on a certain graph: QWs can be used for universal quantum 
computation, this being provable for both the continuous 8| and the discrete version j^. 
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3. 



Experiments have been designed or already performed to implement the QW 

In this paper, we concentrate on the discrete-time QW on a line. We perform a sys- 
tematic study making use of Wigner functions, which are defined for this problem. Wigner 



functions 
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2l| were introduced as an alternative description of quantum states. They 



play an important role in quantum mechanics, having been widely used in quantum optics to 
visualize light states. From the experimental point of view, they provide a way for quantum 
state reconstruction via tomography and inverse Radon transformation I22I . 

Wigner functions are quasi-probability distributions in phase space, meaning that they 
cannot be interpreted as a probability measure in momentum and space configurations. 
This is an obvious fact for any quantum description, and only marginal distributions can 
be associated to probabilities in position or momentum (or any linear combination, i.e. 
any quadrature). In fact, Wigner functions can take negative values, thus invalidating a 
direct link to a probability distribution. This caveat, however, turns out to be a potential 
advantage, for it can be used to identify "true" quantum states. More precisely, the volume 
of the negative part of the Wigner function, its negativity, has been suggested as a figure of 



n.„t to ,u^tify the degtee of Tht. idea h.. been teeently exploited 

|24j | to directly estimating nonclassicality of a state by measuring its distance from the closest 
one with a positive Wigner function. 

When dealing with the discrete QW, one has to account for the extra degree of freedom 
(in addition to the spatial motion): the coin. We consider the simplest case of a two 
level coin. Therefore, the Wigner function has to incorporate this extra index and, with 
the prescription we use, it turns into a matrix. We will propose a rather straightforward 
extended definition of negativity for this Wigner "function". Then the question arises, what 
kind of states of the QW are nonclassical? Does this quantumness increase in time, as the 
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QW evolves through its unitary evolution? We want to explore these questions using the 
Wigner function. A different topic, although it is closely related to the previous one, is 
whether this nonclassicality can be related to the entanglement between the walker and the 
coin, since this quantity is also evolving during the QW evolution. 

This paper is organized as follows. In Sect. II we review the main definitions pertinent to 
the QW on a line. Sect. Ill introduces the Wigner function for our problem and the main 
associated properties. We present some examples showing our numerical results for the 
Wigner function evolution in Sect. IV. In Sect. V we define an extension of the negativity 



to the QW, based on the proposal made in 



23| for a scalar function. We end in Sect. VI 



by summarizing our main results and conclusions. 



II. DISCRETE-TIME QW WALK ON A LATTICE 

The discrete-time QW on a line is defined as the evolution of a one-dimensional quantum 
system following a direction which depends on an additional degree of freedom, the coin 
(or chirality), with two possible states: "left" \L) or "right" \R). The total Hilbert space of 
the system is the tensor product Hs ® He, where Hg is the Hilbert space associated to the 
lattice, and He is the coin Hilbert space. Let us call T_ (T+) the operators in Hg that move 
the walker one site to the left (right), and \L){L\ , \R){R\ the chirality projector operators 
in He- The QW is defined by 

U{e) = T_ ® \L){L\ C{e) + \R){R\ C{e), (1) 

where C{9) = cr^cos^ + (TxSmO, and a^, ax are Pauli matrices acting on He- For 6 = 
71 /A the operator C{6) becomes the Hadamard transformation. The unitary operator U{6) 
transforms the state in one time step as 

m+i)) = umm)- (2) 

The state at time t can be expressed as the spinor 

00 

\m) = E 

n=— o 

where the upper (lower) component is associated to the right (left) chirality, and 
{| n)/?! G Z} is a basis of position states on the lattice. A basis in the whole Hilbert space 
can be constructed as the set of states {| n,a) =\ n)® | a) /n E Z] a = L, R} . 
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a„(t) 
bn{t) 



in) 



(3) 



III. WIGNER FUNCTIONS FOR THE QUANTUM WALK 



An important tool in some fields related to quantum physics is the use of quasi-probability 
distributions. Wigner functions constitute the major example, although other functions as 



the Glauber-Sudarshan P function 
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26| or the Husimi Q function 27| are commonly used 



in quantum optics. For the case of a one-dimensional system with continuous position x and 
conjugate momentum p, the Wigner function is defined as {2^: 



W{x,p) 



ip*{x — y)ip{x + y)e~ 



(4) 



where il){x) is the wave function for the system in a state and we are using units such 
that h= 1. A number of properties can be derived from the definition, the most important 
ones giving the probability in position (momentum) space as the marginal distributions 
obtained by integration over momentum (position), respectively. We refer the interested 



reader to references 
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22| for an overview about these properties. 



We are interested in describing the QW with the help of Wigner functions. The case of 
a finite lattice, with periodic conditions, has been widely studied 
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Here we want 



to describe the QW on an infinite lattice of equally-spaced positions. A proposal for the 
Wigner function to study this problem will be published elsewhere [32|. As discussed in this 



reference, due to the discreteness of the phase space, one needs to double the phase space 
in order to fulfill the necessary properties of the Wigner function. This doubling feature is 
a characteristic of discrete Wigner functions 
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33] . and we will return to this point later. 



Secondly, we need to incorporate the chiral (or spin) degree of freedom. To do this, we 
consider the walker as a spin 1/2 particle moving on the lattice. This analogy allows us 
to make the connection to the extensive literature of Wigner_ functions describing particles 
with spin. They have been extended to relativistic particles [34] and widely used in kinetic 



theory nuclear physics 



or to descri 



3e neutrino propagation in matter 



this approach, wave functions as defined in 



Within 



32| are to be replaced by (Dirac) spinors. Since 



we are interested in a nonrelativistic description, we simply use 2-dimensional spinors and 
define: 



vr ^ — ' 



2ikl 



(5) 
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In the latter equation, \ip{t)) represents the state of the QW at time t. This definition 
can be extended to the case when the state of the system (walker plus coin) is described by 
the density matrix p{t). In such a case, we would have 

W^p{n,k,t) = -e'^-\2{l,a \ p(t) \ n - I, (3)e-^^''K (6) 

71 ^ — ^ 

In what follows, we will omit the chirality subindices, so that W{n, k, t) will represent an 
hermitian matrix in chiral space. Also, it will be understood that summations with natural 
indices are performed over all integers in Z. Substitution of Eq. ([3]) into Eq. ([5]) immediately 
gives 

W(n,k,t) = ie-Ve-- f "'W'"-'**) V (7) 

Notice that the Wigner function is defined over the phase space associated to the Hilbert 
space of the lattice. This phase space is given by pairs {n,k), with ri G Z and k G [— 7r,7r[. 
These coordinates are not to be confused with the position and momentum of the real 
system, although they are associated to them. To avoid ambiguity, we will refer to n and k 
as the phase space spatial and momentum coordinates. 

Using the above definitions, there are many properties that can be proven in a straight- 
forward way. The most important ones are given below. 

By performing the integration over k one readily obtains, for even values of the spatial 
phase space coordinate, 

a„(t) |2 a„(t)6;(t) 



Ama.«)* = 2 I""'*' I ™" . (8) 
J \ 6„(«)<.(i) 1 6„(«) |2 ; 



while, for odd values 



7r 

J W{2n + l,k,t)dk = 0. (9) 

— IT 

This result is a consequence of the above mentioned doubling feature. From here, the 
spatial probability distribution can be recovered by performing the trace over coin variables: 

^Tr I J W{2n,k,t)dk j =| a„(t) p + | b„{t) 1^= P(n,t), (10) 



where P{n,t) stands for the probabihty of detecting the walker at position n, regardless 
of the coin state. 

To obtain the distribution in momentum space we start from Eq. ([3]) and introduce the 
quasi-momentum basis {| k)/k G [— 7r,7r[} (restricted to the first Brillouin zone), which is 
related to the spatial basis {| n) /n G Z} via Fourier discrete transformation, i.e.: 

{n\k) = ^=e*"^ (11) 
V 27r 



After projecting over a given \k) one obtains 



{k I m) = I ''^] I , (12) 

bk{t) 



where 



a,(t) = 4= E e"^"'«n(t); hit) = ^ E e~™'&n(t) (13) 



are the chirality components of the wave function in momentum space. By introducing the 
closure relation for the basis of states {| k)}, one can relate the Wigner function to the 
matrix elements of p{t) in momentum space: 

Wo.p{n,k,t) = - re^"(''-'=)(g,a I p{t) \ 2k~q,/3)dq, (14) 

with I k,a) =1 fc)® I «). Summation over n leads to 

M{k,t) = ^W{n,k,t) =2{k\ p\k), (15) 

n 

where use was made of the equation Yin ^'"'^ ~ 2TT6{k, 2tt), with S{k, 27r) = Ym ^(^ + 2vrm) 
the "Dirac comb" function. Obviously, M{k,t) is also a 2 x 2 matrix. For a pure state we 
have, using Eq. (fT2|) . 



. , aJt) P aJt)bl(t) , 
M{k,t) = 2{k I ^{t)){^{t) I A;) = 2 I ' ^ ^ ' ''^ ' \ (16) 



bk{t)dl{t) \bk{t) 



|2 



The diagonal components of this matrix give (up to a factor 2) the probability in mo- 
mentum space, when chirality is specified, whereas non-diagonal components correspond to 
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coherences between different chiralities. Now, the trace over chirahty provides the probabil- 
ity in momentum space, when chirahty is not measured: 



^Tr{M{k,t))} = ~al{t) + bl{t)^P{k). (17) 

Since the evolution in momentum space is diagonal in k and given by a unitary trans- 
formation (see, for example 38] ), it can be easily proven that the magnitude P{k) defined 
above remains constant with time. 

Another important property that can be derived for the Wigner function of the QW is 
a recursion formula relating W{n, k,t + 1) to other components of this function at time t. 
Using Eq. ([2]) one obtains, after some algebra: 



W{n, k,t + l) = MRWin - 2, k, t)Ml + e-^^'MnWin, k, t)Ml 

+e^'^MLW{n, k, t)Ml + MiWin + 2, k, t)Ml (18) 

where Ml = {\L){L\)C{9) and Ml = {\R){R\)C{9). An immediate consequence of this 
recursion formula is that sites with even n evolve independently of those with odd n. 

IV. NUMERICAL RESULTS 

In what follows, we will discuss a couple of examples showing the main features of the 
Wigner function. We have numerically simulated the QW evolution for various initial states, 
and explicitly computed the Wigner function. We take the lattice large enough, so that 
boundaries do not need to be considered. In practice, this is equivalent to assuming an 
infinite lattice. In the cases we will consider here, the initial state is such that at t = 
the Wigner function is non-vanishing only at phase space points with even n. Then, Eq. 
(fTSj) warrants that W{2s + l,k,t) = 0, Vs G Z, at any time. Therefore, we only plot the 
Wigner function over the part of the phase space with even spatial coordinate. Fig [U shows 
Wjiji{n, k,t) for t = 500, with initially localized conditions, 

\ ij{0)) = ^{\ R) + 1 \ L))® \ 0) (19) 
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Figure 1: (Color online): Wrji component of the Wigner function with initial conditions given as 
in (|19p . after 500 iterations. 



The initial Wigner function for this state can be easily evaluated, with the result 



2n 



I 1 



Since the matrix defining the Wigner function is Hermitian, this component has no imag- 
inary part. The evolved component Wrr after 500 iterations is shown in Fig. [H and the 
time evolution of this component can be seen on Fig. [2l One can observe an intricate 
structure, arising from interference effects. Notice, for example, the similarity with the 



threads mentioned in 



30|. It is interesting to mention that, although the Wigner function 



expands in space, as the walker distribution broadens, it keeps the same structure. The rest 
of components of the Wigner function show a similar appearance. As an example, we have 
represented in Fig. [2] the real part of the off-diagonal component Wrl for t = 500, starting 
from the localized state Eq. (fT9|) . 

The momentum distribution is obtained from Eq. (fT5|) . We show in Fig. |U as an example, 
the RR component of M{k, t) plotted for t = 50, 100, 200 and the same initial condition (fT9|) 
as before. Since the moment k is bounded, the distribution becomes more intricate as the 
QW evolves. At later times, the figure shows more oscillations, although the envelop remains 
constant, in accordance to the self- resemblance of the Wigner function as time increases. 

For comparison, we investigate the Wigner function for a different initial condition. It 
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Figure 3: (Color online): Real part of the Wrl component of the Wigner function with initial 
conditions given as in (|19p . after 500 iterations. 

corresponds to a "Schrodinger cat" in two positions | ± a) which are entangled with two 
chirahties, i.e.: 

\^{0)) = l^{\a,R)+t\-a,L)). (21) 
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Figure 4: MptR^kji) component with the same initial conditions as before, corresponding to times 
t = 50,t = 100 and t = 200 (left to right). 





Figure 5: (Color online): Wrr component of the Wigner function with initial conditions given as 
in (|2ip . and a = 10, after 500 iterations. 

The initial Wigner function for this state is now 



Notice that a = reproduces the locahzed state described above. The evolved component 
Wrr of the Wigner function is shown in Fig. O Fig. reveals the time evolution of this 
component. As compared to the previous case, it shows an even more complicated structure, 
thus suggesting a less classical state. This suggestion will be confirmed in the next section 
by comparing the negativity for both cases. 




(22) 
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Figure 6: (Color online): Contour plots showing the time evolution of the Wrr component of the 
Wigner function starting from the localized state (|2ip with a = 10. From left to right, the sub 
figures correspond to t = 0, t = 100 and t = 500. 



V. NEGATIVITY OF THE WIGNER FUNCTION IN THE QUANTUM WALK 

Consider the Wigner function, as defined in Eq. (j3j). As already mentioned, the fact 
that this function can take negative values implies that it cannot be considered as a classical 
probability distribution. Therefore, non positivity of the Wigner function can be interpreted 
as a measure of the non-classicality of the system. In quantum optics, this is interpreted as 
a signature for non-classical states of light, caused by a quantum interference phenomenon. 
In the context of finite dimensional systems, this idea is exploited in__|39| to establish a 
criterion for entanglement in a system of two spin 1/2 particles. In j40|, a connection is 
found between entanglement and the negativity of the Wigner function for hyperradial s- 
waves. Such states, in D = 2d dimensions, can be interpreted as the wave function of two 
entangled particles in d dimensions. 



A quantitative measure of non-classicality is given by the negativity, as defined in 
In the continuous case, with variables x and p, this volume can be written as: 



23|. 



/CO noo noo noo 

/ [I W{x,p) I -W{x,p)]dpdx = / I W{x,p) I dpdx - 1. (23) 
OO >/ — OO J —OD J —OO 

In deriving the latter equality, we made use of the fact that the total probability is 
normalized to one, so that f^^W{x,p)dpdx = 1. We will use S(W), adopted to the 
discrete case, to compute the negativity of the Wigner function, and to explore whether we 
can relate the deviations from a classical behavior of the QW to the entanglement between 
the walker and the coin. 
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In our case, the Wigner function is a 2x2 matrix, and variables x, p are to be replaced 
by n and as defined in Sect. III. We propose, as a possible generalization of the above 



equation 



32|: 



/n/2 /•7r/2 
[\\Win,k)\\ -Tr{W{n,k))]dk = y2 \\ W{n, k) \\ dk - 1, (24) 

'7r/2 „ J-7r/2 



r.7r/2 /•7r/2 

[I I Win, k) 1 1 -Tr{W{n, k))]dk = V / 

t/2 „ J-k/2 

and where we omitted the dependence on time to simplify the notation. Again, the latter 
equality arises as a consequence of normalization. As a measure of the norm | | of a matrix 
A, we adopt the trace norm, defined as: 

II A 11= Tr^AU, (25) 

where denotes the hermitian conjugate of A. If Ai(n, A;), \2{n, k) are the eigenvalues 
of W{n,k) for a given n and fc, one obtains ||W^|| =| \i{n,k) \ + | X2{n,k) \. In this way, 
Eq. (p^ adopts the following form: 

-77/2 



S{W) 



f y2{\Xi{n,k)\ + \X2{n,k)\-Xi{n,k)-X2{n,k))dk, (26) 
which can be regarded as a natural generalization of (l23|l . 

We have numerically calculated the negativity ( IMl) as a function of time, for the same 
initial conditions considered in the previous section: the localized state and two different 
cat states (Fig. [7]). 

We immediately observe that one obtains a higher degree of "quantumness" for the 
Schrodinger cat states, as compared to the localized state. This was expected from the 
higher degree of complexity observed in the Wigner function for the cat state, and implies 
a larger amount of interference effects. Higher values of the initial separation yield a larger 
negativity. 

The question arises whether this higher degree of quantumness will also imply a higher 
entanglement between the walker and the coin. To this end, we use the entropy of entan- 
glement as a quantity to characterize this property. To be more precise, we compute the 
quantity 

S{t) = -Tr{p,it)\og,pM}, (27) 
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50 100 150 200 250 300 



Figure 7: (Color online): Negativity, as a function of time, for the localized state o = (blue, 
bottom curve) and for the Schrodinger cat state with a = 4 (red, middle curve) and a = 30 (green, 
top curve). 




Figure 8: Entropy of the localized state a = (blue, bottom curve) and for the Schrodinger cat 
state with a = 4 (red, middle curve) and a = 30 (green, top curve). 

where Pc{t) = Trs{p{t)} is the density matrix for the coin, which is obtained by tracing out 
the spatial degrees of freedom. 

Fig. [H] shows S(t) for the locahzed state (a = 0) and for the cat state with different 
values of the (half) separation a. The entanglement entropy is bounded by the (log of 
the) dimension of the coin space: S{t) < 1, and we can see that the time scale to reach 
this maximum is shorter for larger values of a, which also correspond to higher degrees of 
negativity, as seen in previous figures. 
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VI. CONCLUSIONS 



In this paper, we have studied the discrete-time quantum walk on a one- dimensional 
lattice using the Wigner function. We have explored the potential of this phase-space rep- 
resentation to characterize the dynamics of the system. Differently to the case of a scalar 
wave function, we now have a 2 x 2 matrix, defined over a phase space {n,k). We have 
examined two cases, which correspond to different initial conditions. The first case starts 
with an initially localized state at n = 0, as considered in many works in the literature. Its 
Wigner function shows a quite intricate structure, which is build up by interference effects 
on the lattice. We have also considered Schrodinger cat states, in which space and coin 
are initially entangled. It is apparent, from the Wigner function plots, that such states are 
capable of building up higher interference effects, and thus describe states that are even 
"less classical" than the previous one. In order t o q uantify this effect, we have computed 
the negativity associated to the Wigner function |23'|, as defined by the negative volume of 
the function in phase space. This definition has been extended to our case by the use of the 
trace norm of the Wigner matrix. As expected, the cat states give rise to a larger extent of 
negativity, which is even larger when the initial separation of the cat state is increased. In 
accordance to these ideas, one also obtains that the coin-walker entanglement evolves faster 
for the latter states. Altogether, the time evolution of the QW translates into an evolving 
state which separates from classicality. This separation, as time goes on, was expected from 
the fact that the walker distribution expands faster (with a quadratic deviation a ~ t ) than 
its classical counterpart (the random walk, for which cr ~ y/t). 
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